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Abstract 

The method suggested by Balian and Brezin for treating angular momentum reduc- 
tion in the Faddeev equations is shown to be applicable to the relativistic three-body 
problem. 

x. 

1 Introduction 

The Faddeev equations ]T[] provide a formulation of the quantum mechanical three-body prob- 
lem as a compact kernel integral equation. Compact kernel integral equations are well suited 
to numerical solution because the kernel of the equation can be uniformly approximated by 
a finite dimensional matrix, reducing the dynamical equations to a system of linear algebraic 
equations. Numerical solutions of the Faddeev equations require that the abstract operator 
equations be expressed as integral equations formulated in a chosen basis. The equations can 
be simplified by using bases that exploit the rotational and translational symmetries of the 
problem. The method suggested by Balian and Brezin has proved useful in applications. 
This paper illustrates how the Balian-Brezin method generalizes to the relativistic case. 
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2 The Non-Relativistic Case 



The Hamiltonian for a system of three particles interacting with two-body interactions has 
the form 

H = H Q + Y,Va, (1) 

where H is the kinetic energy operator for three particles and = ® h is the two-body 
interaction between particle i and j imbedded in the three particle Hilbert space. In what 
follows the interaction Vu is denoted by V3, and likewise for cyclic permutations of 1,2, and 
3. 

The three-body 2-cluster to 2-cluster transition operators are 

T ij = V j + V l R{z)V j , V 1 = Vj , (2) 

where R(z) '■— (z — H) -1 with z — E ± i0 + . The second resolvent identity implies that the 
transition operators satisfy the coupled integral equations: 

T\z) = V j + J2 V k R k (z)T k \z) , (3) 

where Rk{z) := (z — Hq — V^) -1 is the resolvent of the Hamiltonian for the interacting pair 
plus spectator. 

Equations (|3|) have the same kernel as the Alt-Grassberger-Sandhas equations|3[], which 
are used instead of the Faddeev equations in most numerical applications. 

They are also representative of the type of equation to which the Balian Brezin method 
can be applied. The analysis that follows is valid for any type of connected kernel three-body 
equation. 

In applications Eq. (|3|) is replaced by a linear system of algebraic equations for the matrix 
elements of T^(z) in a chosen basis. The complexity of the algebraic equations is reduced 
by choosing a suitable basis. 

The Hilbert space for the three-body system is the tensor product of three single-particle 
Hilbert spaces. Let pj, rrij, Sj and /ij denote the momentum, mass, spin, and magnetic 
quantum number of the z-th particle. Plane wave basis vectors for the three-body system 
are 

IP1/U1P2/W2P3M3) (4) 

with normalization 

3 

(pi//lP2/J2P3/*3|pi//iP2/4P3/4> = II 6 (Pi ~ Pi )<W^ • ( 5 ) 

%=\ 

The kernel of Eq. ([3]) 

K{%) := VMz) (6) 

commutes with the total linear momentum operator and the linear momentum operator of 
the spectator particle, i. It is separately invariant under rotations of the spectator and the 
interacting pair. 
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The number of non-zero matrix elements of the kernel can be reduced by evaluating it 
in a basis that exploits these symmetries. One basis commonly used is the basis of angular 
momentum eigenstates defined by the following linear combination of the basis vectors (HI), 



\Pqiki; JfijLSjls) 



dqidkilPqiki] ^ i fi j fi k )Y L ^ L (^)Yi m (k i )(sjfi j Skfik\sfis) 
x{lnisii s \jp j ){s i Hijv j \Sns){LnLSns\JHj) , 



(7) 



where 



|Pqikj; fiifizfis) = J dp 1 dp 2 dp 3 \p 1 ^2^3/^3) 5 (pi - Pi(P, qi, k i)) 
x5{p 2 - p 2 (P, q.i, ki))5(p 3 - p 3 (P, q i; ki)) 
and the momenta P, q^, and kj are related to the single particle momenta by 

P = Pi + P2 + P3 , 

k; = ^p, - , 



and 



q/ 



rrii 



m jk 

M Pl M Pjk • 



(8) 

(9) 
(10) 

(11) 



where pj k : = Pj + Pk, M = mi +m 2 + m 3 , rrijk = rrij + m k) and k are cyclic permutations 
of 1, 2, 3. With this choice the Jacobian of the variable change {P, q i; kj} — > {p 1; p 2 , p 3 } is 
unity. 

The matrix elements of the kernel of Eq. ([3]) in this basis are 
(PqA; JtijLSjlslViRiWlP'&h J'n'jL'S'j'l'sT) 



^ (p _p %-C 



SjJ'Sfiju'jSjfSLLfSsS' 



x(kijls\ViRi z 



2M 



M 



2iriim jk 



Wjl's') 



(12) 



where 



is 



Vi 





f p2 






ViRi | 




q 2 M \ 




* 2M 


1 2mim jk J 




p2 


2 M 




-l 


2M 


2rriim jk 


2rrijm k 





(13) 



and Vi is the two-body interaction. 

The basis (0) exploits the translational and rotational invariance of the system and the 
spectator particle. The coupling implicit in Eq. ([3]) breaks the translational and rotational 
invariance associated with the original spectator particle. This can be illustrated by con- 
sidering the iterated kernel. Let K(i) denote the kernel with particle % as spectator and let 
\%) denote the basis vector ([7|) corresponding to particle i being the spectator. The matrix 
elements of the iterated kernel are 



(14) 
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where in order to express the kernel in a basis where it has the form fllzD it is necessary to 
introduce the change of basis Since j ^ i this changes the spectator and breaks the 

invariance associated with the spectator particle. 

The change of basis is computed from the definitions by transforming the bras and kets 
in Eq. (|I2D to the tensor product of one-body basis vectors and evaluating the overlap. The 
general result is determined by taking cyclic permutations of 

(1|2'> = 5(P-P0/dqidk 1 dq 2 dk 2 5(q 1 -q 1 (q 2 ,k 2 ))5(k 1 -k 1 (q 2 ,k 2 )) 
x Y^ L Y^ (kx) (q 2 ) Y lVi (k 2 ) 

x(Jfij\L^ L S^ S )(Sfis\Sl^lj^j)(j^j\lfJ'lSfi s )(sfi s \s2fi2S3fi3) 

x(s3/i3Si/ii|/^)(/VisX|jV^)(s2/i 2 jVj|5V 5 )(LV i 5V 5 |JV J ) , (15) 

where the sums over all repeated magnetic quantum numbers are implicit. This expres- 
sion includes four two-dimensional integrals over the angles associated with each relative 
momenta. Two of the integrals can be done in terms of the angular parts of the delta 
functions. This still leaves two two-dimensional integrals over the angles and two delta func- 
tions that fix the magnitude of the relative momenta. In this form the change of basis is a 
complicated transformation to implement numerically. 

Balian and Brezin exploit the rotational invariance of the overlap matrix elements (|i~5"D 
to facilitate the computation. Rotational invariance implies that the matrix element ( |T5D 
is equal to kronecker deltas in 5jj'5 /JJ/Jj , multiplied by an expression independent of It 
follows that ( |T5j ) is equal to its average over fij. The averaging makes the integrand in ( |T5] ) 
invariant under simultaneous rotations of the vectors qi, ki, q 2 , and k 2 which implies that 
the integrand depends only on the independent invariants qf, k\ and qi • kj. The quantity 
qi • kx is fixed in terms of qf, q%, and k\ or kf, q%, and k\. The result is that, after the delta 
function fixes the integral over the cosine of the angle between qi and ki, the integral over 
the remaining angles is equal to a phase-space factor times the invariant integrand. The 
integrand can be evaluated by choosing the vectors qi, ki, q 2 , and k 2 to have any convenient 
values consistent with the kinematics: 

(1|2') = S(P - P')5(E - E') JJ ' 

87r 2 m 23 mi3 ^ ^ ^ ,^ , 

mim 2 m3/ci/c 2 gig 2 L ^ L iw fJ ' L Mi 

x(Ju\Lfi L SfI S )(Sfl S \SlfJ'ljfij)(jfij\lHlSfi s )(sfi s \s2fi2S3^3) 

X(s3/i3Sl/il|sV s )(ZV/sVjiVi)( s 2/i2//ii|^Vs)(^Vi^Vsl^^ , (16) 



X 



where E is the kinetic energy and the unit vectors are fixed by kj and q^, which can be 
chosen arbitrarily subject to the constraints that fix qf, kf, and kj • qj. Note that the above 
expression includes a sum over the overall magnetic quantum number, v. 

Balian and Brezin suggest three choices of fixing the direction of kj and qj, subject to 
the constraint on the scalar product, which facilitate the computation of these expressions. 
These three choices lead to additional simplifications. This method for treating the change of 
basis matrix elements is employed in many existing numerical treatments of the three-body 
problem. The same considerations apply to configuration space computations. 
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The form of the recoupling coefficient in Eq. (|16|) uses one of the delta functions in the two 
momentum variables to do the integral over kj qj while the other gives the energy-conserving 
delta functions. In principle these delta functions can be used to perform integrals over any 
pair of variables. In applications it is convenient to use the delta functions in the relative 
momenta to do the two integrals over the relative momenta. In this case the one-dimensional 
integral over kj • must be calculated numerically. 

3 The Relativistic Case 

The Balian-Brezin method is also applicable to large class of relativistic formulations of 
the three-body problem, which includes all generalized Bakamjian-Thomas |3J] formulations, 
including specific Bakamjian-Thomas formulations in any of Dirac's|5|] forms of the dynamics. 

Relativistic invariance of a quantum model requires that all probabilities have values 
independent of the choice of inertial coordinate system, where by definition any two inertial 
coordinate systems are related by Poincare transformations continuously connected to the 
identity. Wigner's theorem states that this condition implies the existence Q of a unitary 
ray representation of the Poincare group on the three-body Hilbert space. 

The problem of constructing the unitary representation of the Poincare group plays 
the same role in relativistic quantum mechanics as constructing the unitary representation 
of the one-parameter time evolution group in non-relativistic quantum mechanics. In the 
non-relativistic case it is sufficient to solve the eigenvalue problem for the Hamiltonian. 
The one-parameter time evolution group is expressed in terms of the eigenfunctions and 
eigenvalues of the Hamiltonian as 

U(t) = £ \E)e- iEt (E\ . (17) 

E 

In the relativistic case it is sufficient to solve the eigenvalue problem for the Casimir operators 
of the Poincare group. For system of massive particles this is equivalent to finding the 
simultaneous eigenstates of the mass and spin operators. Representation theory for the 
Poincare group is used to construct the unitary representation, U(A,a), of the Poincare 
group in terms of the eigenvalues and eigenstates of the mass and spin similar to the manner 
that the eigenvalues and eigenstates of H are used to construct U (t) in Eq. ( |17D . 

Bakamjian-Thomas models are models with interactions that commute with and are 
independent of three independent functions of the non-interacting four momentum and 
all components of a chosen non-interacting spin. The analysis that follows is limited to 
Bakamjian-Thomas type models. The dynamical problem in a Bakamjian-Thomas dynam- 
ics is to find simultaneous eigenstates of the mass and non-interacting spin operator in a 
suitable basis. The Balian-Brezin method can be applied to Bakamjian-Thomas models be- 
cause the mass operator commutes with the non-interacting spin vector, however, the explicit 
justification in the relativistic case is more complicated because the quantities that replace 
qi and kj in the relativistic case do not generally transform as vectors under rotations. 

The mass operator in the relativistic case has the same form as the Hamiltonian in the 
non-relativistic case, consisting of a mass operator M for three non-interacting particles 
plus pairwise interactions 

M = Mo + £>y- (18) 

i<j 
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Time dependent scattering theory |7| gives the transition operators in the relativistic case, 

T ij := V J + V*R{z)V j , V 1 = Y,Vj , (19) 

where R(z) = (z — M)~ l is the resolvent of the mass operator. The second resolvent identity 
implies the components of the transition operators satisfy the coupled equations 

T*(z) =Vi + J2 V k R k (z)T k ^z), (20) 

where for the relativistic equation Rk{z) := (z — Mq — \4) _1 is the resolvent of the mass 
operator for the interacting pair plus spectator. Eq. fl2"0|) has the same form as the non- 
relativistic equation. Differences appear in the relation of the two-body interactions in the 
three-particle Hilbert space and the two-body interactions in the two-body Hilbert space. 
The structure of the interactions is given in Eq. fl43f) below. 

The Hilbert space for a single particle in relativistic quantum mechanics is an irreducible 
representation space for the Poincare group corresponding to the mass and spin of the parti- 
cle. Vectors in this space can be expanded as linear combinations of simultaneous eigenstates 
of the linear momentum and ^-component of spin. In relativistic quantum mechanics there 
are many different spin operators that satisfy SU(2) commutation relations whose square 
is the total spin. Although the spectrum of the magnetic quantum numbers is fixed by 
the SU(2) commutation relations, the physical interpretation of the spin operator is de- 
termined by the transformation properties of the single particle states under the Poincare 
group. In general the single particle vectors |p/i) transform as mass-m spin-s irreducible 
representations of the Poincare group 

1/2 

D s ^(B-\A Pm )AB(p m )) , (21) 

where a> m (p) := VP 2 + m2 is the energy, p m = (o> m (p), p) is the four momentum of a particle 
with mass m and momentum p, Pa := Ap m is the Lorentz transform of the four momentum 
p m , and B{p m ) is a Lorentz-boost-valued function of p m with the property 

B(p m )p = p m , (22) 

where po := (m, 0, 0, 0) is the 0- momentum four- vector. The quantity 

^m(PA) 
W m (p) 

fixes the normalization of the transformed state to ensure the unitarity of U(A, a). 

The combination R w (A,p m ) := B^ 1 (Ap m )AB(p m ) is a rotation for any A, called the 
Wigner rotation associated with the boost B(p m ). The interpretation of the magnetic quan- 
tum number is determined by the choice of B(p m ). A boost B(p m ) is defined up to a rotation 
valued function R(p m ) of p m 

B(p m ) -> B\p m ) = B(p m )R(p m ) . (24) 



t/(A,a)| P/U ) = e iAp -- a |p A /i / ) 



Wm(p) 





<9pA 




dp 



(23) 
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The canonical boost B c (p m ) is the unique boost with the properties 

BM = I (25) 

and 

R w (R,Pm) = B' 1 (Rpm) RB C (p m ) = R (26) 

for any rotation R. Eq. (p6|) states that for canonical boosts the Wigner rotation of a 
rotation is the rotation itself. This is not true for other types of boosts. 

Any other boost is related to a canonical boost by a p-dependent rotation as in Eq. (p4|). 
This rotation is called a generalized Melosh rotation^]. The helicity spin and front-form 
spin are examples of spins corresponding to non-canonical boosts|7j]. These two choices are 
distinguished by other special properties. 

A basis for the three-particle Hilbert space is the tensor product of three single particle 
bases. Basis vectors have the form 



|PlyUlP2Ai2P3A i 3) (27) 

with the same normalization convention as the non-relativistic expression (|]). 

A basis for the three-body system that simplifies the matrix elements of the kernel K(i) = 
ViRi(z) of the relativistic three-body equations (^Cj) is constructed by finding the coefficients 
of the linear transformation that take the product of three irreducible representations of the 
Poincare group to a direct integral of irreducible representations. This is equivalent to the 
problem of constructing Clebsch-Gordan coefficients of the Poincare group. This replaces 
the successive pairwise coupling of the irreducible representation spaces of the Euclidean 
group (rotations and translations) used in the non-relativistic reduction. 

The construction of the Clebsch-Gordan coefficients of the Poincare group is most easily 
understood by successive coupling of pairs of irreducible representations. To treat the general 
case first consider the case of canonical spin. Kinematic variables for the two particle system 
are 

P = Pl +p 2 , m 12 = ^[^P~P , ki = B-\P) Pi , (28) 

corresponding to the total four momentum of the non-interacting pair, the invariant mass of 
the non-interacting pair, and the relative momentum of the non-interacting pair, respectively. 
The relative momentum defined in Eq. (|28|) is not a true four-vector. It undergoes Wigner 
rotations when the system is Lorentz-transformed. 

The two-body basis vectors that transform irreducibly under the action of the tensor 
product of one-body representations of the Poincare group are 

\Pfi(kjls)) := Z7i[S c (P)]®^ 2 [S c (P)]|dk|k,/ii,-k,/i 2 )^ w (k) 
x (s 1 fi 1 s 2 n 2 1 sp, s ) {luisn, \jy) 

^milPlJ^m 2 lP2j rn 12 



J «ik|p 1 (P,k),/i' 1 ,p 2 (P,k), y u 2 ) 



1/2 



^mi(ki) UJ m2 (k 2 ) to 

mi2 V x J 

wc (B c (P),k 2 )]Y l n(k) 
x(si^is 2 ^ 2 |s/x s )(//i Z s/i s |j» , (29) 



xD^R^B^M^DlJRUBmMW^i 
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where k = ki = — k 2 . The factor m 12 /c<j mi2 (P) fixes the normalization to be 

if Mk-k>) 



(p/x(fcj/s)|py(fcyz's')> = s(p - p' 



k 2 



Sjj'5uu,'5u'S ss i . 



(30) 



The single particle and combined spins can be transformed from canonical to any other 
type of spin (i.e., a spin associated with an arbitrary boost B x {p m )) with the unitary 
transformation |7j 

\ Pf i) x = |pz/)^yiUp)] , RcM ■= B- 1 (p)B x (p) , (31) 
where R cx (p) is a generalized Melosh rotation. A direct calculation using the relations 

(32) 

and 



Y lm (k c )D l mi/l [R xc (p)] = Y lvi (R xc (p)k c ) = Y ln (K) 
R xc {pi)R wc {B c {p), k lc )R 

cxip) — Rwx\Bx\P)i k\x) Rxc(M-lx) 



2) 



(33) 



shows 

\Pl*(kjls)) a 



dk|pi(P, k), n[, p 2 (P, k), (jf 2 ) 3 



UmAPl) U m2 {p 2 ) 77112 



(k 2 ) 



1/2 



xD^JR wx (B x (P), k lx )R xc (k lx )}D%jR wx (B x (P), k 2x )R xc (k 2x )} 



xy iw (ki a; )(si//iS2A<2|s^ s )(^sA*sb» 



(34) 



In expression ( |3"4] ) the argument of the D-functions is the product of a generalized Melosh 
rotation followed by the Wigner rotation associated with the x-spin. The quantity k^ is 
related to the canonical k c by a generalized Melosh rotation 

k x = B~ 1 {p)B c {p)k c = R xc (p)k c . (35) 

As a result of this construction the state \Pfi(kjls)) x transforms like a particle with mass 
rriij = \J mf + k 2 + J m 2 + k 2 and spin j where the magnetic quantum number transforms 
as an x-spin: 



Ui(A, a) ® Uj(A, a)\PfjL(kjls)) x 

W m,,( P A) 1 " 



JAP-ai 



PKjj!{kjls)) 3 



^my (P) 



^[^(A,P)]. 



(36) 



The transformation properties of this two-particle state are identical to those of a single 
particle. 

To construct the three particle basis the tensor product for a pair plus spectator 



\PjkHjk{kdls)) x <g> |pi//i) 5 



(37) 



is decomposed into a direct integral of irreducible representations by repeating the above 
analysis replacing one of the single particle states by the two-body state ([34]). The new 
kinematic quantities are 



P = Pi+ Pj 



M n 



V^P~P , q l = B- 1 (P) Pi , q jk = B- 1 (P)p jk , (38) 



8 



where q := q« = — q,fc- In addition, the two-body kinematic variables are 



kj = B c 1 (p jk )p j , h = B c 1 (p jk )p k , k := k, = -k* 
The resulting basis is related to the single particle bases by 



(39) 



|P/i; qJLSkjls) 



x 



dqdk|pi^Pj/4p fc /4) 

VmjiPj) Um k (Pk) m jk 



1/2 



Umjfe) w mfe (k fc ) W mjfc (Pjfc) u mi (q_i) u^iqjk) ^Afo(P) 
X-D*/ w [R W x(B x (p jk ), kj^Rxdkj^D 8 ^ [R wx (B x {p jk ), k kx )R xc (k kx )} 

x y, w (kk,) (sjfijSkHk I s// fl ) (Z^s// S | j 2/) 

X -^ > l'.Uj [Rwx(B x (P), qj kx )R xc (qjkx)]Dp w [P^P^P), fe)-Ra;c(fe)] 

xY L ^ L (cL ix )(jVj Si/iilS ^(L^lS/isI J Vj) ■ (40) 

This basis also transforms irreducibly under the tensor product of three one-body represen- 
tation of the Poincare group, 



C/i(A,a) <g) U 2 (A,a) ® U 3 (A,a)\P/i(qJLSkjls)) x 



e lP ^ a \P^\qJLSkjls)) x 



^Mo(P) 



P^[P_(A,P)] 



(41) 



The unitary representation of the Poincare group in Eq. ([|T]) is not the physical representa- 
tion. What is relevant is that in a Bakamjian-Thomas formulation both the interacting and 
non-interacting representation of the Poincare group share the same spin operator. 

There is a collection of Bakamjian-Thomas models associated with each different type of 
irreducible basis element of the form ( 4"0|) that are unitarily equivalent and have the same 
scattering matrix elements. The kernel of the integral equation is defined in terms of its 
matrix elements in the appropriate basis 



(Pp, qJLSkjls\ViRi(z)\P'n'; q' J'L'S'k'j'l's 1 ) 



5(P - P^Sjjrf 



5(q - q') 



huhs^AkMV^R^z-q^k'jl's') 



(42) 



where the operators Vi(q 2 )Ri(z; q 2 ) are related to two-body interactions Vi with matrix ele- 
ments (jkls\vi\jk'l's') by 



q- + (mjk + Vi) 

2 



q- 



m 



'jk 



q- 



q 2 + (m jk + Vi) 2 



(43) 



As in the non-relativistic case, this can be computed only after solving the two-body problem. 
What is important is the collection of delta functions on the right hand side of Eq. (0). For 
Bakamjian-Thomas models, where the functions of the four-momentum that commute with 
the interaction are not the three components of the linear momentum, the delta functions 
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in the components of the linear momentum must be replaced by delta functions in the 
appropriate functions of the non-interacting four-momentum. For instance, in Dirac's point- 
form dynamics |5|] the three components of the four- velocity replace the three components 
of the momentum. In the case of Dirac's front-form dynamics^ the delta functions in the 
components of the three momentum are replaced by delta functions in the three components 
of the four momentum that generate translations tangent to the light front, x 3 + t = 0. The 
front-form case also requires a special choice of spin operator. 

As in the non-relativistic case, when Eq. (|20|) is iterated, a new interaction is introduced, 
violating the symmetries associated with the initial spectator particle. Thus if (i\K(i)\i) 
denotes the expression in Eq. (|42| ) the matrix elements of the iterated kernel are 

In order to take advantage of the simple form of the kernel in Eq. ( f4"2| ) it is necessary to 
compute the overlap where 

(1|2'> := (P//; qx JLxSxhjihsxlP'v'; q' 2 J> L' 2 S' 2 k' 2 f 2 l' 2 s' 2 ) , (45) 

and cyclic permutations. By direct calculation the matrix element is 

3 

<1|2'> = /rfqA /dq / 2 (ik / 2 ^5[p ^ (P,q 1 ,k 1 )-p i (P^q^,k' 2 )] 

i=l 

y u m2 (p 2 ) ^ m3 (p 3 ) m 23 <j mi (pi) u m23 (p 23 ) M 

UJ m 2 (k 2;c ) ^m 3 (k 3;c ) UJ m23 (P23) W mi (qiz) ^m 2 3 (<123:e) &M (P) 

1 /2 

^m 3 (P 3 ) ^mi(p'i) rn' u m m2 (p' 2 ) u m ' 3i (p' 31 ) Mq 

(p'31) (q'si X 

x(J / u|L/i L S , /i 5 )(S'/i5|j/i i Si/ii)y L * Mi (qi a; ) 

xDJ HUj [R cx (q 23x )R wx (B-\P),p 23 )]D^ 
x(7'i/j|Z^js/i a )(s^ 8 |s2/x 2 s 3 /i3)y { * | (ki a .) 

xD T z v< 3 [ R wx{ B x{'P'zi)i Kx) R xc{k' 3x )}Dl\^ [R wx (B x (p' 31 ), k' lx )R xc (k' lx )] 
xY i'^ ( k 2 J (ssv'zSiis'xls' )i' s ) (I'^s'^lj'v'j) 

xD^, [R WX (B X (P'), q' 3lx )Rxc{q', lx )]DH^ [R WX (B X (P>), q' 2x ) R xc (q' 2x )} 
><^V L J (q 2 J0'V>2^|5y 5 )(LK^|jy) . (46) 

Following Balian and BrezinJJ, symmetry principles are used to evaluate this matrix ele- 
ment. To facilitate the evaluation of the matrix element (1|2') note that 

(P/x; qxJLtSthjihstlP'ii'; q' 2 J' L' 2 S 2 k 2 j 2 l' 2 s 2 ) 

= (P/i;g 1 JL 1 S 1 A; 1 j 1 / 1 s 1 |^ t (A,a)f/(A,a)|Py;g 2 J , L 2 5 2 A; 2 j 2 / 2 s 2 ) (47) 

for any A and a. The kinematic quantities are 

ki = k 2x . = -k^, k' 2 = k' 3x = -k' lx , qix = -q 23a: , c( 2x = -q[ 31x (48) 
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with 



k 2x = B x 1 {p 23 )P2, hx = B x \p 23 )p 3 , k' 3x = B x \p' 3l )p[ 
and 

q lx = B; 1 (P)p 1 , q 23x = B x -\P)p 23 , q' 2x = B x \P)p' 2 , 
Evaluating the right hand side of Eq. (|47|) gives 

wm„(Pa)^(P^ 12 



3; 



k' lx = B^ipM (49) 
Qsi x = B- l (P)p' 3V (50) 



(1|2'> 



JA(P'-P)-a 



^Mo(P) ^M'(P') 



D J V ;[R WIB (A,P))D J V ,JR^^,P')] 



x (P A u; q 1 JL 1 S 1 k 1 j 1 l 1 s 1 \P' A p'; q' 2 J> ' L' 2 S 2 k f 2 f 2 l' 2 s' 2 ) . 



(51) 



To deduce the implications of Eq. (|51~1) consider different choices of A and a. For A = I 
and arbitrary a the relations (|51| ) cannot be satisfied unless P = P' . For A = B X (P), 
P = P' = P = (M , 0, 0, 0), and a = it follows that 



(P/i; giJIriS'ifciji/iSilPV; q' 2 J' L' 2 S' 2 k' 2 j 2 l' 2 s' 2 ) 

1/2 

(P /i; gi JLi^i^uiZiSilP'o/i'; q' 2 J' L' 2 S' 2 k' 2 i' 2 l' 2 s' 2 ) (52) 



M 



M 



o; Mo (P) ^(PO 

which shows that Eq. ([51]) can be expressed as 8(P— P') multiplied by a quantity independent 
of P. For the case that A = R is a rotation, P = Pq = RPq, and a = it follows that 



(P /x; gi JLxSifciji/iSilPo//; q 2 f L 2 S 2 k 2 j' 2 l' 2 s 2 ) 
= (P i/; gi JLiS'ifcijiZiSi|Pdi/; q' 2 f L 2 S 2 k' 2 j' 2 l 2 s 2 ) 



x B'l*[R wx (R, Po)]D^, ,[R U1X (R, P )} . 



(53) 



The left hand side of Eq. (|53|) is independent of R. Since the rest boosts B X (P ) are rotations 
(normally chosen to be the identity), R wx (R,Po) is a representation of SU{2) which can 
be parameterized by elements of 577(2) and integrated over the group. Since the Haar 
measure P| is normalized to unity the integral is still equal to the left side of Eq. (|53|). The 
integral over the D functions can be done explicitly using M 



dRD u *^[R wx (R, Pq)]D u ,^ [R wx (R } P )\ 
with the result that 



1 



2J + 1 



(54) 



(P /i; qiJLiS^xjihs^PofM; J' q 2 L' 2 S 2 k 2 j' 2 l' 2 s' 2 ) 
J J 2J+1 



(P i/; giJLi^i^ijiZiSilP'o^; J' q' 2 L' 2 S' 2 k' 2 i' 2 l' 2 s 2 ) 



(55) 



The general form of the matrix element is obtained by incorporating all of the conse- 
quences of Eq. (pl|): 



(P/i; g x JLiSifcijiZiSilP /i ; g 2 J L 2 S' 2 k 2 j 2 l 2 s 2 ) 
= 5jj,5^6(P - P / )5[M (q 1 , kO - M (q' 2 , k 2 )] 
xA J (giLiS'i/ciji/iSi;g 2 L 2 S' 2 /c 2 j 2 / 2 4) 



(56) 
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where the amplitude Aj(l;2) is invariant under the non-interacting representation of the 
Poincare group. From the definition (|56f) the amplitude Aj(l; 2) can be computed from the 
matrix element (1|2') in Eq. ([46]), 



8{M - M')Aj(l; 2') 

= I J ^PiPiiiqxJLiSiktjthsilPo^q^JL^ky'^s^) . (57) 

Since Aj(l; 2') is invariant, the computation of this quantity is facilitated by evaluating this 
expression for P' = Pq. 

The last step used by Balian and Brezin is to exploit the rotational invariance. To do 
this in the relativistic case let the total four momentum be Pq = RPq and evaluate 

(P /x; qiJLxStkijihs^P'oii] q' 2 J' L' 2 S' 2 k' 2 j' 2 l' 2 s' 2 ) 

= J (Po/i;giJL 1 5 , i/CijViSi|^ t (i?)f/(i?)|pi/iip 2 /i2P3^3)rf 3 Plrf 3 P2^ 3 P3 
x(p 1 /i 1 P2^2P3^3|^ t (i?)^(i?)|P^;g2^2^^2J2^4) 

= J D^*(R)(P u]q 1 JL 1 S 1 k 1 j 1 l 1 s 1 \Rp 1 /j, 1 Rp 2 lJ,2Rp3^3)d 3 pid 3 p 2 d 3 p 3 
x (Rp m Rp 2 ^Rp m \Py; q' 2 J'L' 2 S' 2 k' 2 f 2 l' 2 s' 2 )D J u ^(R) 

(P v; giJL 1 S , i/cij'i/iSi| J Rp 1 /ii J Rp 2 /i2 J Rp3/U3)^ 3 Pi^ 3 P2C? 3 P3 
x {Rp^Rp^Rp^P'^- q' 2 f L' 2 S' 2 k' 2 i' 2 l' 2 s' 2 ) . (58) 

This equation means that the matrix element (1|2') averaged over the magnetic quantum 
numbers is invariant under simultaneous rotations of all of the single particle momenta that 
appear as intermediate states. The rotated qj's and kj's obtained by rotating the single 
particle momenta are (for P = P ) 

qra = Rwx(R, Po)<ii , k ffi = R WX (R, gjfc)ki . (59) 



To apply the above result to the computation of Eq. (f4~6D, note that there are nine 
delta functions and eight variables of integration. Three of the delta functions lead to the 
overall momentum conserving delta function, leaving six delta functions. Since the k do 
not transform simply under rotations, it is practical to use four of the remaining six delta 
functions to perform the integrals over ki and k 2 . Two delta functions remain. One factors 
out of the expression, giving the conservation of the invariant mass. The other can be used 
to fix the angle between qx and q 2 . Three integrals over the two unit vectors qi and q 2 with 
fixed qi • q 2 remain. If the matrix element is averaged over the magnetic quantum numbers 
then the resulting integrand is necessarily independent of the remaining three variables of 
integration. The integral is 8n 2 multiplied by the integrand. The result, after computing all 
of the Jacobians needed to convert the delta functions to the desired form, is: 



(1\2') = S(P-P')5(M-M r S "~ 1 



■ 

2 



X 



2J+1 |k la .||ky|q la .||q 2a 



Um 2 (kix Vm 3 (ki x )o; mi (qi x )a; m2 3 (qi x )a; m 3 (k 2x )u mi {k! 2x )uj m2 (q!2 X ) UJ m 31 (q' 2x 



1/2 
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x ( J/j, I Lfi L Sfi s ) ( Sfi s | j n j a i \i i ) Yl^ L ( q i x ) 
xDl. v .[R^(q 23x )R wx (B-\P ),q 23x )]D^ 

X (Pj I SUslfll) (SH S I S 2 /i 2 S 3 /i 3 ) (kia,) 

x£>^, [R WX {B X {P£), QsiJRMiJWU [R wx (B x (Po), <4Rc(gL)] 

x>LV L J (q 2 J0'V>2/i 2 |^)(LV^V 5 l^) • (60) 

The invariant part of Eq. (^) has been evaluated with P = P , which implies the replace- 
ment of all of the p^s by the corresponding q^s. In the case of the front form or point form 
the spin must be chosen accordingly (front-from spin for front-form interactions, canonical 
spin for the point-form interactions) and the delta functions in Eq. (p0|) are replaced by 

8(P - P')S(M - M') -> 5(P + - P +/ )5(Pi - Pi)S(P 2 - P 2 )5(M - M') (61) 

in the front form, and by 

5{P - P')5(M — M') — > 5{y — V) ^ M ~ M ') ( 62 ) 

in the point form, where {P + , Pi, P 2 } are the front-form components of the four-momentum 
and V = P/M are the independent components of the four- velocity. 

Note that if the magnitude of all of the k's and q's appearing in both the rotation 
matrices and the Jacobian factor in Eq. ( |B"0"| ) are set equal to zero this expression has the 
non-relativistic expression (|16|) as a limit provided the rest boosts are chosen to be the 
identity. Specifically the D functions all become the identity and the factor inside the 
| • • • I 1 / 2 becomes mi3m 2 3/m 1 m 2 m 3 . The resulting expression is identical with Eq. (|T6"D. 

The expression ( |60D explicitly involves the four unit vectors qi, q' 2 , k 1; and k' 2 . The 
q- can be evaluated in any geometry, subject to the constraint that the angle between the 
two unit vectors is fixed by kinematic considerations. These choices fix the quantities k,. 
Different choices of the geometry used to evaluate the q- can lead to additional simplifications 
in the evaluation of Eq. (|60D, although in the relativistic case the choice of best geometry 
depends to some extent on the choice of spin. If the boosts satisfy B X (P ) = I then the rest 
Wigner rotations R^B' 1 (P )q)] in Eq. fl60D can be replaced by the identity. For canonical 
spin there are no generalized Melosh rotations, while for the front-from spin the Wigner 
rotations of front-form boosts are the identity, R w f(Bf{p),q) = I. All of these properties 
lead to further simplifications of Eq. (|60|). 

Each of the choices of spin and continuous variables in a Bakamjian- Thomas model im- 
plies a choice of representation. Although it is tempting to formulate a model using a choice 
that leads to the simplest Racah coefficient, any choice has implications for the structure 
of the representation of other operators, such as electromagnetic and weak current opera- 
tors. The interacting system has interactions in different Poincare generators for different 
choices of the Racah coefficients. For instance, the Racah coefficients leading to an instant- 
form dynamics imply that the infinitesimal generators of rotationless Lorentz transformation 
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contain interactions, while Racah coefficients appropriate for a point-form dynamics imply 
interactions in the momentum operators. For more general choices of Racah coefficient there 
can be interactions in any number of Poincare generators. For an operator, such as an elec- 
tromagnetic current operator, that is well approximated in an impulse approximation in one 
representation may require large two-body contributions in another representation. Thus, 
other considerations may be important in choosing a representation. 

As in the non-relativistic case, it is useful to replace the delta function in the kinematic 
mass by a delta function that expresses the invariant mass explicitly in terms of the relative 
momenta, 

S( M - M') - S(M(k u qi ) - M(k' 2 , q> 2 )) , (63) 

where in computing this delta function it is important to note that only three of the variable 
can be considered independent. This replacement allow one to do the integral over one of 
the relative momenta. 

The analysis above shows that the methods suggested by Balian and Brezin is appli- 
cable to relativistic three-body equations of the Bakamjian-Thomas type. The relativistic 
expressions for the recoupling coefficients were shown reduce to the non-relativistic ones in 
the non-relativistic limit. 
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